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Abstract. Spatially-resolved gas pressure maps of the Coma galaxy cluster are obtained from a mosaic of XMM- 
Newton observations in the scale range between a resolution of 20kpc and an extent of 2.8 Mpc. A Fourier 
analysis of the data reveals the presence of a scale-invariant pressure fluctuation spectrum in the range between 
40 and 90kpc and is found to be well described by a projected Kolmogorov/Oboukhov-type turbulence spectrum. 
Deprojection and integration of the spectrum yields the lower limit of ~ 10 percent of the total intracluster 
medium pressure in turbulent form. The results also provide observational constraints on the viscosity of the gas. 
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1. Introduction 

In hierarchical structure formation scenarios clusters grow 
via accretion and merging of smaller subclumps. Gas ac- 
creting onto clusters of galaxies has bulk velocities of 
about v = 1900(T/6.7keV) a52 kms- 1 at 1 Mpc (e.g. 
Miniati et al. 2000), where T is the mean X-ray tem- 
perature of the intracluster medium (ICM). This veloc- 
ity is comparable to the expected sound speed of 1000- 
1500 km/s of the ICM. Accretion flows through filaments 
and sheets are highly asymmetric and produce complex 
patterns which can survive for long time-scales in the ICM 
(Miniati et al. 2000). Simulations by Norman & Bryan 
(1999) predict that the turbulent pressure in the ICM can 
account for up to 20% of the thermal pressure. We thus 
expect some measurable effects of turbulence in the ICM 
of clusters of galaxies. 

Concerning X-ray data, Inogamov & Sunyaev (2003) 
propose a study of spectral line profiles as a useful di- 
agnostic tool of turbulent flows in the ICM which could 
be measured with the future ASTRO-E2 satellite. Vogt & 
EnBlin (2003) propose the application of Faraday Rotation 
measures to test turbulence in the ICM, and claim that 
for a few clusters a Kolmogorov spectrum seems to be 
plausible. 

In the present investigation we show that turbulence 
in the ICM can be probed directly with pressure maps 



Send offprint requests to: Peter Schuecker 
peters@mpe.mpg.de 

* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by 
ESA Member States and the USA (NASA) 



provided by the XMM-Newton satellite as a result of its 
high sensitivity and excellent spectral capabilities. 

Section[21 summarizes the basic phenomena related to 
turbulent flows. In Sect.[5]we give a simple analytic treat- 
ment of projection effects introduced through observation. 
In Sect.^we present the observational data and describe 
how our X-ray temperature and pressure maps are con- 
structed. Based on the direct comparsion of local temper- 
ature and density measurements, we give in Sect. |3] some 
arguments that their observed fluctuations appear to be 
almost adiabatic. The same statistical analysis also sug- 
gests the absence of pronounced contact discontinuities 
and strong shocks. These observations provide a baseline 
consistent with the presence of a turbulent flow. Therefore, 
we study in detail in Sect. El the measured pressure spec- 
trum in Fourier-transformed k space and discuss its inter- 
pretation in Sect.0 For all computations a flat geometry 
and a Hubble constant of H — 50 km s _1 Mpc -1 are used. 
We assume a distance of 139 Mpc to the Coma cluster so 
that 1 arcmin corresponds to about 40kpc. 

2. Phenomenology of isotropic turbulence 

Traditionally, the phenomenology of isotropic turbu- 
lence is based either on second-order velocity statistics 
(Kolmogorov 1941) or on their Fourier-transformed coun- 
terparts (Oboukhov 1941). The velocity energy spectrum 

C K t 2 ' 3 k- 5 / 3 , 



E v {k) 



(1) 



with the non-dimensional Kolmogorov constant Ck, can 
be obtained for the inertial scale range from simple dimen- 
sional arguments (e.g., Lesieur 1997). In Eq. @J, E v (k) 
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Fig. 1. Angular profile of the weight function (/3 profile 
squared) integrated along the z direction for various integration 
lengths. The lowest/highest contours indicate 10/100% contri- 
butions, and the intermediate lines the contributions in steps 
of 10%. The upper panel shows the contributions on cluster 
scales and the lower panel the contributions on the maximum 
scales of the observed turbulence structures (see Sect.|7J. For 
the brightness distribution we assume a f3 profile with the pa- 
rameter values f3 = 0.75 and r c = 420 kpc. 
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Fig. 2. Pressure power spectra with an intrinsic slope of 
n = — 7/3 as expected for a Kolmogorov/Oboukhov turbulence 
(r c = 0), and its projection (420 kpc), as seen along the z di- 
rection through the central plane of a cluster with a core radius 
(r c = 420 kpc) and slope parameter (f3 = 0.75) as measured for 
the Coma cluster. In order to illustrate the projection effects 
over a large scale range, we did not introduce any characteristic 
scale which limits the spectra at large and small scales. 

Tab. 1. Spectrum E v (k) in the inertial range for different flu- 
ids, Mach numbers M, and magnetic fields B. 
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B-field 
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v{k) 
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< 1 
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-5/3 
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> 1 


B = 


k 


-6/3 
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< 1 


B^0 
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-5/3 


[4] 


> 1 


B^0 


k 


-3/2. ..-9/3 



[1] Kolmogorov (1941), Oboukhov 1941), [2] Burgers (1974), 
[3] Goldreich & Sridhar (1995), [4] Cho & Lazarian (2002), 
Vestuto et al. (2003). 



is the kinetic energy of fluctuations per unit mass and 
wavenumber k (physical units m 3 s~ 2 ), and e, in units 
of kinetic energy per unit mass and time (m 2 s -3 ), is 
the rate of kinetic energy transport from large to small 
scales. Eq. Q fundamentally stems from the assumption 
that, based on dimensional arguments, a turbulent eddy 
on a scale A ~ k , decays on a 'turnover' timescale 
r = X/v(k) = [k^E^k)]- 1 / 2 (Oboukhov 1941). Since the 
kinetic energy associated with fluctuations over a scale 
fc -1 is kE v (k), in steady state regime the rate of energy 
transport across different scales is e = kE v /r which leads 
to Eq. (JTJ. 

In the inertial scale range where Eq. applies, turbu- 
lence develops without being affected by boundaries, ex- 
ternal forces, or viscosity. Here, the fluctuating quantities 
are assumed to be statistically invariant under translation 
(homogeneity) and rotation (isotropy) . Tab. 1 summarizes 
some of the results from theoretical studies and numerical 
simulations, which suggest that Kolmogorov/Oboukhov- 



like spectra emerge in an inertial scale range under quite 
general conditions. 

While all these studies are based on the analysis of ve- 
locity fluctuations, Oboukhov (1949) and Batchelor (1951) 
showed that gas pressure fluctuations also obey a scaling 
law (e.g., Lesieur 1997, Chap. VI), 

E P {k) = C P e^k- 7 l\ (2) 

where Cp is a non-dimensional constant, and Ep(k) has 
the units of kinetic energy per unit mass squared and is 
normalized to unit wavenumber k (physical units m 5 s~ 4 ). 
We note that the slope of the spectrum of the pressure 
is steeper than the spectrum of the velocity (P ~ v 2 ). 
As for the velocity spectrum, we also expect that in gen- 
eral the exact slope of the pressure spectrum depends on 
whether or not the fluid is supersonic, and whether or not 
magnetic fields are present. The pressure spectrum thus 
appears as an excellent and powerful diagnostic tool of 
turbulent ICM flows. In addition, this approach appears 
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quite attractive because pressure fluctuations can already 
be measured from high resolution X-ray data. 

Establishing the presence of a turbulent ICM implies 
testing whether or not there is a scale-free inertial range 
in the pressure spectrum with a slope similar to the 
Kolmogorov/Oboukhov case. It is of further interest to 
measure the location of certain characteristic scales, such 
as the spectral break which corresponds to the scale where 
the kinetic energy is initially injected into the ICM, as well 
as the smallest scale where the corresponding energy is fi- 
nally dissipated into the ICM. It is, however, not yet clear 
whether there is also energy dissipation within the inertial 
range caused by the development of randomly distributed 
weak shocks (Burgers turbulence, Tab. 1). 

3. Projection effects 

X-ray observations measure, after filtering and normaliza- 
tion (Sect.EJ, a projection of the actual three-dimensional 
pressure fluctuations, 6P/P(r), on the two-dimensional 
celestial sphere. Note that the application of normalized 
instead of absolute quantities modifies the physical units 
of the structure function and power spectrum introduced 
in Sect. |21 Our analysis concentrates on scales which are 
small compared to the cluster core radius r c . As larger 
scales of order r c are approached, the global cluster pro- 
file starts to be probed. This will be taken into account in 
the treatment (see Sect-EJ- Thus, we start by decomposing 
the fluctuations into Fourier modes, 



SP 

7W 



— ik-' 



(3) 



The small-angle approximation allows us to regard the 
cluster region as composed of coplanar layers, with a 
pressure distribution assumed to be homogeneous and 
isotropic within each layer. We deliberately choose an X- 
ray energy band which is almost free of the temperature 
dependence of emissivity (see Sect.0J. A two-dimensional 
pressure field can thus be constructed using r = (R, z), 



A 5P ( ^ 

dz — (r) 
P y ' 



1 



(z - zo) 7 



-3/3 



(4) 



which weights the pressure field of each layer with its 
emissivity, as obtained from a f3 model of the cluster gas 
density (squared). In Eq. zq is the distance between 
the observer and the cluster centre along the line-of-sight 
(LOS), and R a two-dimensional vector on the sky. 

An illustration of the weighting scheme is shown for 
the Coma cluster in Fig-QJ Here, the contours give the 
percentage of surface brightness contributed by ICM gas 
within a distance \z — zq | from the plane through the clus- 
ter centre and perpendicular to the z direction, as ob- 
tained by the integral of the squared fj profile along the z 
direction. A comparison of the upper and lower panels of 
Fig. fl] reveals that the angular dependence of the general 
profile of the Coma cluster is imprinted on large scales, of 
order > 200 kpc, whereas on smaller scales the structure 



of Coma appears quite homogeneous in each layer, so that 
only a unique profile along the z direction is seen. In the 
inner region, weighting can thus be approximated by tak- 
ing into account only the variation of the density along 
the z direction and neglecting the angular dependence as 
assumed in Eq. . For clusters with small core radii the 
approximation is less valid and the projection should be 
performed numerically. 

A direct consequence of Eq. (TIJ or more complicated 
projection schemes is the invariance of relations between 
fluctuating quantities under geometric projections. This 
is illustrated for the adiabatic relation between tempera- 



ture and density, T ■ 



1 , used in Sect .0 to classify the 



fluctuations. Its differential version 

(r) = 7 - 1 — (r . 

1 n 



(5) 



defines the (adiabatic) fluctuations. The projected tem- 
perature and density fluctuations are given by Eq. J2J re- 
placing P with T and n. We further replace [1 + (z — 
z o) 2 / r c]~ by W(r), i.e., a general weighting function 
which describes the geometric projection process. If the 
weighting W(r) were the same for T and n, then we could 
write 



/ST 
dz — (r) W(r) 



f(H, 



= (7— 1) const / dz — (r)W(r) 
J n 

= (7 — 1) — (R) ■ 

n 

However, two-dimensional (projected) temperature maps, 
T(R), resulting from X-ray observations, are n 2 (r)- 
weighted averages along the LOS of the three-dimensional 
field T(r). On the contrary, projected squared densities, 
n 2 (R), are obtained by a simple geometric mean without 
any weighting. For this more realistic case we thus have 



ST(R) _ J dzST(r)n 2 (r) 7 - 1/ dzT(r)Sn 2 (r) 



T(R) ~ J dzT(r)n 2 {r) 



J dzT{r)n 2 {r) 



,(6) 



where 5n 2 = 2n Sn and Eq. JSJ) have been used. For the 
assumed central core region of the cluster we now have 
the advantageous situation that T(r) assumes the role of 
an almost constant weighting function for the density, and 
we can show that retaining only first order terms, Eq. JSJ 
simplifies to 



ST(R) 7-1/ dz5n 2 (r) 7-I Sn 2 (R) 



T(R) 



J dzn 2 (r) 



n 2 (R) 



(7) 



Eq. 10 is the observational counterpart of Eq. (J5J and will 
be used in Sect .[S] to constrain possible types of pressure 
fluctuations measured in the centre of the Coma cluster. 

After cross-correlating projected fluctuation measures, 
we now proceed with projecting pressure spectra. We still 
assume that the central cluster region, to which our anal- 
yses are restricted, is small enough to be approximated 
by a set of coplanar, homogeneous, and isotropic layers, 
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Fig. 3. Projected temperature map (upper left), pressure map (upper right), entropy map (lower left) and image substructure 
(residual) map as revealed by small-scales in the wavelet decomposition (lower right). The maps are obtained from spectral 
hardness ratios and surface brightness data and smoothed with a wavelet filter of the Coma cluster. Each map covers an area 
of 93 x 93arcmin 2 . 



yielding a unique pressure profile along the z direction. 
We want to take advantage of the invariance property of a 
Gaussian profile under Fourier transformation by replac- 
ing the f3 model by a Gaussian with variance r 2 /3/3. The 
approximation is better than 5% for r < r c and allows us 
to regard the convolution of the pressure field with the gas 
density profile in Eq. 10} along the z direction as a trans- 
formation 5P/P{k) -> 5P/P(k)exp(-k 2 r 2 /60). With 
the standard relations between two and three-dimensional 
power spectra (e.g., Peacock 1999, Sect. 18.1) we obtain 
the simple expression 



fluctuations along the z direction leads to an exponentially 
damped two-dimensional spectrum. Damping is larger for 
smaller modes because relatively more fluctuations are 
added along the z direction. 

The relations between power spectra and pressure 
spectra in two and three-dimensions can be obtained from 
the condition that the sum over both statistics must give 
the same total energy, 



1 



1 



dk z 7> 3D y/kl + K 2 ) exp 



3/3 



(2irf 
1 



2nKdKV 2 T>(K) = [ 
Jo 

/ 4irk 2 dkT 3 B(k) = / E 3B (k)dk, 
Jo Jo 



E 2D {K)dK, (9) 



where the integration extends over the wavenumber k z 
along the z direction. The wavenumbers in two and in 
three dimensions are K and k, respectively. Note that 
K = 2tt/Q, with being the angular scale of the pro- 
jected fluctuation, has the physical units rad - , but can 
be readily transformed to metric scales using the distance 
to the cluster. Eq. (jHJ shows that adding three-dimensional 



(8)0)3 

which yields the following conversions: 

V 2 b{K) =2t<K- 1 E 2 d(K), 
V w {k) = 2n 2 k- 2 E 3D (k). 

This gives the projection of the pressure spectrum 
dk z 



E 2D (K)=K 



k 2 +K 2 



E m Wkl + K 2 ), 



(10) 



(11) 
(12) 



(13) 
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Fig. 4. Detailed view of the projected pressure distribution of the central region of the Coma cluster. The 145 kpc scale 
corresponds to the largest size of the turbulent eddies indicated by the pressure spectrum (SectJTJ. The smallest turbulent 
eddies have scales of around 20 kpc. On smaller scales the number of photons used for the spectral analysis is too low for reliable 
pressure measurements. 




Fig. 5. Nested grids of temperature (left panel), pressure (middle panel), and entropy (right panel) measurements. Each map 
covers an area of 69.3 x 69.3 arcmin 2 . 



6 



Schuecker et al.: Probing Turbulence 



where 2?3d( - ) can be identified with the spectrum Eq. (J2J 
after including appropriate upper and lower integration 
limits. Note that the decline of (k) and thus of Vsd (k) 
outside the inertial range may not be sharp, and instead 
better described by smooth characteristic shape functions 
(see, e.g., Fritsch et al. 1980). Fig. [21 gives a quantitative 
impression of the projection effects on a power spectrum 
computed with Eq. J5J for the Coma cluster with distance 
139 Mpc, j3 = 0.75, and core radius r c — 420 kpc. For illus- 
tration, no K limits of the inertial range are introduced. 
For comparison we also plot the k z r c <C 1 case where 
no corrections for projection are necessary. The intrinsic 
slope of the pressure spectrum is n = — 7/3. Note that the 
power spectrum of the pressure fluctuations is damped at 
60 kpc by a factor of about 28. The observation of pressure 
fluctuations along the z direction through the cluster thus 
damps their amplitudes by a factor \/28 ~ 5.3. However, 
cluster-wide fluctuations on Mpc-scales would appear al- 
most undamped. 

Further observational effects are related to the mea- 
surement errors of the pressure and the global pressure 
profile and can be illustrated best with the observed spec- 
trum of the cluster pressure distribution (Sect.HJ). 

4. The Coma mosaic 

In this paper we use the performance verification observa- 
tions of the Coma cluster obtained with the EPIC-pn in- 
strument on board of XMM-Newton (Jansen et al. 2001). 
Previous reports of these observations were given by Bricl 
et al. (2001), Arnaud et al. (2001), Neumann et al. (2001, 
2003) and Finoguenov et al. (2004a, for point-like sources). 
This work includes all datasets obtained to date, as de- 
scribed in detail by Finoguenov et al. (2004a, see also 
Tab. 2). While most of the pointings have been obtained 
in Extended Full Frame Mode (Frametime 199 msec in 
Tab. 2), three observations were conducted in Full Frame 
Mode (Frametime 73 msec). 

All observations have been reprocessed using the lat- 
est version of the XMM reduction pipeline (XMMSAS 
5.4.1), which yields an astrometry to better than 1 arcsec. 
Although the Coma data are public, for some (Coma-10 
and Coma-0) of the early observations of the performance 
verification phase no complete Observational Data Files 
(ODF) had been produced by standard processing, and a 
special preprocessing (XMMSAS task odffix) was done on 
such pn exposures at MPE by Michael Freyberg. As a re- 
sult, a few Coma pointings are not yet publicly available, 
which precluded us from using MOS data. The vignetting 
correction, crucial for obtaining reliable source character- 
istics over a wide region, is performed using the latest cal- 
ibration (Lumb et al. 2003). Two pointings at the Coma 
centre were used in that calibration, in a way requiring 
that the same sky pixels yield the same flux between the 
two observations. The level of the emission was not used 
in the calibration, so it could be analyzed further. The 
RMS fluctuations of the comparison of two Coma fields is 
within 2%, which will affect the apparent pressure fluctu- 



ations studied here on the 1% level, much lower than the 
observed 10% amplitudes. 

The images were extracted separately for each point- 
ing, along with the corresponding exposure maps. 
We select pn events with PATTERN < 5 and 
(FLAG&0xc3b0809) = 0, which in addition to FLAG = 
events includes events in the rows close to gaps and bad 
pixels; however, it excludes the columns with offset energy. 
This event selection results in a better spatial coverage of 
the cluster, but at a somewhat compromised energy reso- 
lution, which is sufficient for the broad-band imaging. 

Our final results are derived from a spectral analy- 
sis where only the FLAG = events were retained. For 
background subtraction we used the similarly screened 
and selected events from the background accumulation of 
Andrew Read (Read & Ponman 2003) and also subtracted 
out-of-time events as a background, using products from 
the SAS task epchain. This subtraction is important as 
some pointings (see FrameTime 73 ms in Tab. 2) are per- 
formed in the Full Frame Mode. 

To provide an overview of the structure of the ICM 
of the Coma cluster, we show in Figs.|3 and 0] the tem- 
perature, pressure, and entropy maps, as well as maps of 
the small-scale surface brightness structure. These maps 
use hardness ratios in the 0.8-2 keV and 2-7.5 keV bands, 
calibrated for the measurement of temperature, as a sub- 
stitute for the temperature determined directly from the 
spectral analysis. The projected entropy (S) and projected 
pressure (P) maps are derived from the projected tem- 
perature, T, and surface brightness, E, through the re- 
lations S = TIT 1 / 3 and P = TE 1 / 2 , respectively. The 
maps are constructed from composite wavelet filtered im- 
ages to suppress the large scale background. The details 
of the analysis based on surface brightness and hardness 
ratio maps and the rationale of the use of wavelet filtering 
is described in detail in previous publications (Briel et al. 
2003, Finoguenov et al. 2004b, and Henry et al. 2004). 

The substructures seen in these maps suggest 
turbulent-like fluctuations. The pressure maps are of spe- 
cial importance because they clearly show fluctuations 
which are not contaminated by contact discontinuities (see 
Sect.EIl ■ For a quantitative study of the significance of 
these fluctuations we thus performed direct fits to the 
spectral X-ray data. 

In Fig.[5]we show the temperature, entropy and pres- 
sure maps, based on the temperature and emission mea- 
sure obtained through direct spectral fitting, using sev- 
eral grids to define the region of spectral extraction. 
Only the 16 x 16 grids with a pixel size 40 x 40 arcsec 2 , 
120 x 120 arcsec 2 , 260 x 260 arcsec 2 are shown. This figure 
also illustrates the relative positioning of the grids. Each 
binning involves a mixing of various spectral components. 
Therefore, a decision has to be made on which of the com- 
ponents the spectral analyses should be performed. We 
have chosen to put our interest on the hotter component, 
and so have used the 1-7.9 keV energy band for spectral 
fitting. Fine grids, with a pixel size of 40 x 40 arcsec 2 and 
lower, located in the central region, do not suffer that 
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Name of pn-Camera observing times (ksec) Orbit Frametime 



Date Observation RA(2000) DEC(2000) planned effective msec 



2000 May 29 


CoiTlcl C6Ilt6r 


12 


59 


46.7 


27 57 


00 


15.0 


12.8 


86 


199 


2000 June 21 


CoiTlcL 1 


12 


56 


47.7 


27 24 


07 


25.0 


21.1 


98 


73 


2000 Tune 1 1 


Coma 2 


12 


57 


42.5 


27 43 


38 


25.0 


43.8 


93 


199 


2000 Tune 27 


W 111 CI • 1 


12 


58 


32.2 


27 24 


12 


25.0 


11.0 


101 


73 


2000 June 23 


Conici 4 


13 


Of) 


04.6 


27 31 


24 


25.0 


5.4 


99 


199 


2000 Mav 29 


nnmfl Pi 

V_. Wllltl fj 


12 


59 


27.5 


27 46 


53 


20.0 


9.9 


86 


199 


2000 Tune 12 


V,, W 111 CX vj 


12 


58 


50.0 


27 58 


52 


20.0 


12.4 


93 


199 


2000 Dec 10 


Coma 7 


12 


57 


27.7 


28 08 


41 


25.0 


18.6 


184 


199 


2000 Dec 10/11 


Coma 8 


13 


01 


25.6 


27 43 


53 


26.0 


14.3 


184 


199 


2000 June 11/12 


Coma 9 


13 


00 


32.7 


27 56 


59 


20.0 


14.7 


93 


199 


2000 June 22 


Coma 10 


12 


59 


38.4 


28 07 


40 


20.0 


15.4 


98 


199 


2000 June 24 


Coma 11 


12 


58 


36.5 


28 23 


56 


25.0 


11.5 


99 


199 


2002 June 5/6 


Coma 12 


13 


01 


50.2 


28 09 


28 


25.0 


9.6 


456 


199 


2002 June 7/8 


Coma 13 


13 


00 


36.5 


28 25 


15 


25.0 


20.0 


457 


199 


2000 June 22 


Coma (bkg) 


13 


01 


37.0 


27 19 


52 


30.0 


12.8 


98 


199 


2001 Dec 4/5 


Coma cal 


12 


59 


46.6 


27 57 


00 


25.0 


17.4 


364 


73 



much from temperature mixing, but they do suffer from 
small number statistics. So, for those we used the 0.5-7.9 
keV band. A detailed check has shown that for a similar 
location in the Coma cluster all grids yield similar tem- 
perature estimates, which supports our choice of energy 
bands. The selection of the grid resolution was performed 
to yield at least 5000 counts per pixel. The total number 
of counts available for the analysis in the Coma observa- 
tion reaches two million counts in the 0.5-2 keV band and 
a similar amount in the harder band (2-7.9 keV). 

5. General character of the fluctuations 

In order to obtain more information about the type of fluc- 
tuations seen in Figs.OHSl we performed a pixel-by-pixel 
cross-comparison of temperature and density gradients. 

FigurcEl shows the correlation between the gradients 
of the projected X-ray temperature and the gradients of 
the projected squared gas density as obtained for the 40 x 
40 arcsec 2 pixel grid. We concentrate on this specific grid 
because it is mainly restricted to the core region of the 
Coma cluster and has sufficiently high signal-to-noise X- 
ray spectra at a comparatively small angular resolution. 

The relative fluctuations are determined for each pixel 
of the temperature and density map by averaging the gra- 
dients over its four nearest neighbour pixels. The averaged 
gradients are obtained for density and temperature maps 
and can be compared in a point-wise manner. However, 
the individual fluctuations are large (10% level) because 
they include also the measurement errors. Therefore, an 
additional binning with dn 2 /n 2 = 0.1 and a continuous 
sliding of this bin along the density axis is necessary to 
see a clear trend. The error bars are the la fluctuations 
of the mean obtained for each bin. 

For the classification of the fluctuations, we show in 
Fig.H3 model expectations obtained with Eq. (JZJ). The line 
labeled 7 = 5/3 corresponds to a monoatomic ideal gas. 




-0,04 0,04 

<dT/T> 

Fig. 6. Correlation between relative fluctuations of density 
squared (n 2 ) and temperature (T), and their la errors. The 
two thick lines represent the adiabatic exponent 7 = 5/3 
(monoatomic ideal gas) and 7 = 4/3 which gives a good repre- 
sentation of the data. Due to the sliding window method used, 
neighbouring data points are correlated. 

For contact discontinuities local pressure equilibrium leads 
to dn/n = —ST/T. Apparently, the 7 = 4/3 line gives a 
better representation than what is expected for the 7 = 
5/3 case. This could be due to a contamination by contact 
discontinuities. 

Since the gradients are measured on scales of 27kpc, 
which are small compared to the cluster core radius of 
r c = 420 kpc, we expect them to be sensitive probes of lo- 
cal substructure fluctuations and not significantly affected 
by the global cluster profile. We nevertheless tested this 
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approximately with Monte Carlo simulations where the 
gradients of n 2 are computed with an isothermal f3 model 
of the Coma cluster and a gas adiabatic equation of state. 
The gradients are determined in the same way as the em- 
pirical data and added to the adiabatic density fluctua- 
tions. We found that large-scale gradients in the density 
field broaden the scaling relation, but without introducing 
a bias in the determination of 7. In fact, we have verified 
this for different values of the adiabatic exponent 7. 

The observed temperature and density maps are also 
tested for possible correlations between T and n 2 intro- 
duced by the large-scale distribution of the ICM. The 
radial profiles are obtained by averaging temperatures 
and densities in concentric rings with a width of 50 kpc. 
Whereas n 2 shows a significant decrease of 27 ± 2 percent 
between the cluster center and 300 kpc, the temperature 
decreases by only 3± 2 percent relative to the central value 
of 7.34 ± 0.13 keV. The observed temperature gradients 
thus appear with the same size as the errors and can thus 
be neglected as a possible second-order effect. Within this 
approximation, no correlations between temperature and 
density fluctuations are introduced by the global cluster 
profile. 

To conclude, Fig.El suggests a positive correlation be- 
tween temperature and density gradients which is not re- 
lated to the large-scale distribution of the ICM. The gradi- 
ents occupy different regions than contact discontinuities 
and strong shocks. The data are actually quite close to 
the expected adiabatic case. In order to find out whether 
or not such fluctuations are organized as in a turbulent 
regime, we study in the following the statistics of the spa- 
tial pressure fluctuations. 

6. Power spectrum of spatial pressure fluctuations 

The first step in our (standard) power spectrum analysis is 
the determination of the global pressure profile P(R) from 
the observed 2-dimensional pressure map P(R), in order 
to get the residual local pressure fluctuations, 5P/P(R) = 
P(R) / P(R) — 1. The second step is the determination of 
the Fourier power spectrum of SP/ P(R), corrected for the 
errors of the pressure measurements (shot-noise subtrac- 
tion), and normalized to unit number of Fourier modes 
and to unit area in K space. The resulting projected spec- 
trum V2T>(K) has the physical units kpc 2 . In the follow- 
ing example, the pressure is measured in a regular grid of 
32 x 32 cells, each with 20 x 20arcsec 2 . This grid covers 
the central core region of Coma up to 431 kpc and has the 
fundamental mode K — 2n/X = 0.0146 kpc -1 . The results 
obtained with the other three grids are given at the end 
of this section. 

The global pressure profile P(R) is obtained from a 
low-passband Fourier-filter applied to P(R) with a filter 
scale of 150 kpc, which leaves the global cluster profile 
above this scale almost unchanged. To illustrate the effect 
of the filter, we show in Fig.|7|the power spectrum (marked 
'P') obtained from a direct Fourier-transformation of 
P{R)- On scales between 20 and 40 kpc, the spectrum has 
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Fig. 7. Projected power spectra of different angular pressure 
distributions from the 20 x 20arcsec 2 grid. Lower thin con- 
tinuous line: raw spectrum including shot-noise, substructure, 
and cluster profile (P). Upper thin continuous line: spectrum 
from the normalized pressure distribution (dP/P). Thick con- 
tinuous line: spectrum of the global cluster pressure profile 
as determined with the wavelet transform. Dashed horizontal 
lines: shot-noise levels computed from measured pressure er- 
rors. Dashed vertical lines: characteristic scales. 

a flat plateau-like distribution which is determined by the 
temperature and density errors (shot-noise, see below). 
Between 40 and 125 kpc the spectrum increases signifi- 
cantly above the shot-noise level. This is the spectrum of 
the substructures seen in Figs.|3to[SJ Beyond 125-150 kpc, 
the spectrum abruptly increases due to the global pressure 
profile of the Coma cluster. 

A similar increase is also seen in the spectrum marked 
Wavelet' which is obtained alternatively from a wavelet- 
filtered pressure map. For the wavelet decomposition we 
used the algorithm of Vikhlinin et al. (1998) and com- 
puted the spectrum from the wavelet reconstruction of the 
30 x 30 arcsec 2 map with the lowest angular resolution. The 
wavelet algorithm performs a self-adjusting noise suppres- 
sion so that almost no significant shot-noise occurs in the 
spectrum of the global pressure profile. The 20 x 20 arcsec 2 
grid does not cover the complete cluster area and is thus 
not optimal for the proper sampling of the global cluster 
profile. Therefore, the similarity of the 'P' and 'Wavelet' 
spectra is not very good on large scales. However, grids 
with larger bin sizes cover larger scales and give a very 
good agreement with the 'Wavelet' profile (see below). For 
the following analyses we thus use the Fourier low-pass fil- 
ter with a filter scale of 150 kpc to determine P(R) for all 
four grids. 

The histogram of the resulting SP/P(R) is shown in 
Fig-IHl Their distribution appears quite consistent with a 
Gaussian random field (KS-probability of 90%) with zero 
mean and a standard deviation of 15 percent (including 
shot-noise) on a pixel scale of 13.5 kpc. The SP/P(R) field 
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Fig. 8. Comparison of a Gaussian profile with the histogram 
of projected pressure contrasts SP/P = [P(R)/P] — f in con- 
figuration space for a pixel size of 20 x 20 arcsec 2 which cor- 
responds to 13.5 x 13.5 kpc 2 . The mean over all fluctuation is 
almost zero and the la standard variation a$ P / P = 0.15. 



can thus completely be summarized by a power spectrum. 
The corresponding power spectral densities are marked 
by 'dP/P' in Fig.[7| The spectrum shows the expected 
drop beyond 150 kpc. A similar drop at scales below 20 kpc 
marks the resolution limit as given by the pixelation (see 
'Point Sources'). These two cutoff scales limit the range 
of the power spectrum of the substructures. 

The effect of temperature and density measurement 
errors is seen in the power spectra as an almost scale- 
independent shot-noise level which must be subtracted 
from the 'dP/P'-spectrum (dashed horizontal lines in 
Fig-d) ■ For the determination of the shot-noise, we first 
determine at each grid point the local relative pressure 
error op, using T and its error ot as well as n 2 and its 
error a n i , as obtained from the local X-ray spectral fit. We 
then draw at each grid point a random value for the local 
relative pressure error from a Gaussian distribution with 
zero mean and standard deviation a p. After performing 
this randomization for all grid points, one realisation of a 
random map is generated with fluctuations solely caused 
by measurement errors. For the determination of the shot- 
noise we averaged the power spectra of 100 random reali- 
sations. 

The 1 a error bars shown in Fig. are also determined 
from the variances of the spectra obtained from random- 
ized maps of the measurement errors. The errors are lower 
limits because they are obtained from unstructured pres- 
sure maps. Unfortunately, much larger effort is needed to 
improve these estimates, for example, with a set of hydro- 
dynamical cluster simulations. 

Figurc|ni shows the power spectrum of 6P/P(R) after 
shot-noise subtraction. This spectrum can be compared 




0.04 



0.06 



0.08 0.1 
K [1/kpc] 



0.2 



Fig. 9. Projected shot-noise subtracted power spectral densi- 
ties, P2D (K) , of the pressure fluctuations, and their la errors, 
obtained from the 20 x 20 arcsec 2 grid. 



via Eq. (JSJ with theoretical 3-dimensional power spectra 
or structure functions. 

The same analysis performed on the 40x40 arcsec 2 grid 
gives the power spectrum shown in Fig. llQI The errors are 
smaller compared to the results obtained with the 20 x 
20 arcsec 2 grid. We attribute this to the higher signal-to- 
noise X-ray spectra obtained with the larger pixels. The 
spectral shape appears somewhat more curved and steeper 
than the spectrum obtained with the smaller grid. The 
120 x 120 arcsec 2 and 260 x 260 arcsec 2 grids mainly sample 
the global cluster profile. The resulting 'P'-spectra shown 
in Fig.^2 follow the profile obtained with the wavelet- 
filtered pressure distribution. These spectra do not have 
any significant fluctuation power on scales below 100 kpc 
so that we do not show the corresponding 'dP/P'-spectra. 

The differences seen in the power spectral densities 
shown in Figs.lTJl and ITUI are caused by the chosen center, 
bin size, and total size of the sample grid. This sample 
variance can be reduced by averaging the spectral densities 
measured at the same K- values. The average is meaning- 
ful because all spectral densities - although determined in 
K bins with different sizes (fundamental modes) - are nor- 
malized to the same unit volume of 1 /kpc 2 in K space. The 
random errors of the spectral densities which do not refer 
to sample variance are not reduced because the two power 
spectra cannot be regarded as completely statistically in- 
dependent. We regard the spectrum shown in Fig.^1 as 
the final result of the power spectrum analysis. 

7. Discussion 

The present investigation aims to detect turbulence in the 
ICM of the Coma cluster using the pressure spectrum. 
Under certain approximations, one also expects a scale- 
invariant spectrum of temperature fluctuations to be a 
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ties, T , 2d{K), as in Fig.|§lfor the 40 x 40arcsec 2 grid. 
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Fig. 11. Projected power spectral densities as in Fig.0for the 

120 x 260 arcsec 2 and the 260 x 120 arcsec 2 grid. The spectra 
basically follow the global cluster pressure profile. 



probe of ICM turbulence (e.g. Lesieur 1997). However, 
this relies on the assumption that temperature behaves as 
a passive scalar. Once this has been verified empirically, 
the almost uniform distribution of the temperature over 
scales 3> r c allows a cleaner distinction between small- 
scale turbulent substructures and the large-scale cluster 
profile. However, in reality temperature maps are affected 
by cold fronts and other contact discontinuities which con- 
taminate the diagnostic maps. On the other hand, pressure 
maps have a clear relation to velocity (P = pv 2 ) and are 
not significantly contaminated by contact discontinuities. 
Therefore, we regard pressure as a more direct probe of 
ICM turbulence. 
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Fig. 12. Observed projected shot-noise-subtracted power 
spectral densities (dots with la error bars) as obtained for 
the 20 x 20 arcsec 2 and 40 x 40 arcsec 2 grids, compared with 
model predictions (dashed lines). 

The mosaic of XMM-Newton observations is well- 
suited for the detection of turbulence in the Coma cluster 
because it allows a better geometric discrimination be- 
tween pressure variations originating from the overall clus- 
ter profile, and substructure superimposed onto it. This 
transition occurs at about 150 kpc. 

The measured temperature and density gradients 
(Fig.[5J) suggest that the substructures have an adiabatic 
exponent of 7 k 4/3, which is close to the adiabatic case 
of an ideal monoatomic gas. On the other hand, contact 
discontinuities and strong shocks seem to be less likely in 
the core region, consistent with hydrodynamical simula- 
tions (Miniati et al. 2000, Miniati 2003). In addition, the 
statistics of the residual pressure fluctuations appear quite 
Gaussian (Fig.|SJ) emphasising their random nature. Their 
Fourier power spectrum thus completely summarizes the 
fluctuating pressure field and can be used to obtain ob- 
servational evidence for the presence of turbulent flows 
which are characterized by a Kolmogorov/Oboukhov-like 
spectrum. 

Figur eTT2l shows the combined power spectrum of the 
Coma cluster on scales between 40 and 90 kpc. A scale- 
invariant range of the spectrum is indicated and suggests 
the detection of an inertial range of a turbulent ICM. 
Theoretical three-dimensional power spectra, 



V 3 u{k) = 2tt 2 k- 2 E m {k) 
= 2n 2 Ck n ~ 2 , 



27T 2 C P e 4 / 3 k n ~ 



(14) 
(15) 



are transformed with Eq. JHJ) into their two-dimensional 
counterparts (dashed lines in Fig. 1 121) . From the compar- 
ison of observed and theoretical spectra we see that on 
scales between 40 and 60 kpc, the observed power spec- 
trum has a slope between n — —7/3 and —5/3. This 
slope corresponds to the spectral slope of the Fourier- 
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Fig. 13. Contribution of the turbulent pressure to the to- 
tal thermal pressure (contour lines of equal percentage) for 
a Kolmogorov/Oboukhov spectrum with the slope n = — 7/3 
and the amplitude C = C P e 4/3 = 0.0063 kpc~ 4/3 . The power 
spectrum is integrated between the injection and the dissipa- 
tion scale. 



transformed Kolmogorov/Oboukhov structure function 
(Eq.EJl . On scales between 60 and 90 kpc the spectrum 
bends towards smaller slopes between n = —5/3 and 
— 1/3. The normalization constants range from C — 
C P e 4 / 3 = 0.0063 kpc~ 4/3 for n = -7/3, to C = 
0.470 kpc 2 / 3 for n = —1/3. An absolute calibration of the 
pressure used in our analyses is in preparation and will 
give a direct estimate of e. In this respect it would be in- 
teresting to compare this rate of kinetic energy transport 
with the observed X-ray luminosity. 

The power spectrum shown in Fig.^|allows a first esti- 
mation of the location of the characteristic scale where the 
spectrum sharply drops towards larger scales. This scale 
is at approximately Ai rs 100 kpc and should be regarded 
as a lower limit because of possible contaminations by the 
global cluster profile. This scale should also roughly cor- 
respond to the injection scale (e.g., Lesieur 1997), and it 
is similar to estimates for the impact parameter of merg- 
ing clusters based on kinematics and tidal torque-based 
arguments (e.g., Sarazin 2002). 

The integral of the power spectrum (Ea.llO|) is ex- 
pected to give important information about the energy 
deposited in turbulent motion. For the scale range be- 
tween 40 and 90 kpc, the slope and amplitude parameters 
derived above yield relative contributions of the turbulent 
pressure to the thermal pressure between 7.4 percent for 
n = —7/3, and 6.6 percent for n = —1/3. The largest pos- 
sible contributions arc obtained with the n — —7/3 spec- 
trum. Therefore, we computed for this spectrum the rela- 
tive contribution for different minimum scales, i.e., lower 
integration limits of the inertial range (dissipation scale) , 
and maximum scales, i.e., upper integration limits (injec- 
tion scale). 

FigurelT^l shows that for a fixed turbulence spectrum 
the relative contribution is mainly determined by the value 



of the injection scale Ai. We do not see any turbulent ed- 
dies of the size of the core radius of 420 pc in Fig. 0] which 
could have erroneously been subtracted by the Fourier 
low-pass filter - although they still might be present, but 
are difficult to discriminate from the global cluster pro- 
file. Therefore, the relative contribution of the turbulent 
pressure to the thermal pressure should be smaller than 
25 percent. If we take the indication for a turnover in the 
power spectrum shown in Fig.EH at A; = 100 kpc as the 
injection scale, then we would get a lower limit of about 
10 percent. The simulations of Norman & Bryan (1999) 
suggest additional support by turbulent pressure of about 
20 percent, averaged over the cluster (5 to 35 percent be- 
tween core and virial radius), which is apparently of the 
same order as the present observational limit. However, 
further study is definitely required in order to establish 
how the observational quantities relate to the simulation 
results. 

For the observed turbulent ICM we can now esti- 
mate the kinematic viscosity by assuming that magnetic 
fields have a negligible effect (see below). For a turbu- 
lent flow the Reynolds number of the global fluid Sft mea- 
sured at the injection scale Ai, and the Reynolds number 
Sftd measured at the dissipation scale Ad are related by 
SK/SRd = (Ai/Ad) 4 / 3 . The power spectra do not show any 
tendency to decrease at A = 30 kpc (Fig-EJ. Therefore, Ad 
is smaller than 30 kpc. In the following we assume a fidu- 
cial value of Ad = 10 kpc. The turbulent flow in the central 
region of Coma can thus be characterized by 5ft/5ftd > 20. 
Reynolds numbers at dissipation scales are expected to be 
above unity so that 5ft will have values in excess of 20. 
Although this estimate is rather conservative, this is the 
best that can be obtained by direct observations at the 
moment. 

For the viscosity we further need the velocity at the in- 
jection scale. This number can be obtained from hydrody- 
namical simulations (Miniati et al., in preparation), which 
typically give for an 8keV cluster a dispersion turbulent 
velocity of v\ i = 250 km s -1 on scales of A; as 100 kpc. 
This provides a quite reliable upper limit to the kinematic 
viscosity of 



v < 3-10 z 



( V -h ) 
V250kms- 1 / 



A; 



100 kpc 



-1 r 



cm 



.(16) 



Note that the coherence lengths of magnetic fields in the 
cores of galaxy clusters as obtained from Faraday Rotation 
measurements are about 5-10 kpc (e.g., Taylor & Perley 
1993) and thus below the scale range covered by the 
present data. Therefore, we regard the upper limit (Eci. ll6f) 
as not significantly affected by magnetic fields. 

Fabian et al. (2003, see also Reynolds et al. 2004) as- 
sume laminar flow of the ICM around the radio galaxy 
NGC 1275 with v x = 700 kms -1 in the centre of the 
Perseus cluster on A = 14 kpc scales. From the laminar 
appearance of the filaments they assume that the effec- 
tive Reynolds number is less than 1000 so that they es- 
timate the lower limit v > 4 • 10 27 [ £J j-]. The upper limit 



12 



Schuecker et al.: Probing Turbulence 



obtained from a turbulent regime and the lower limit ob- 
tained from a laminar regime can be used to estimate the 
range 10-30 kpc where the transition from a turbulent to 
a laminar flow could occur. This corresponds to a dissi- 
pation scale of the ICM in the same range. A remark of 
caution is, however, necessary here, because we compare 
two different situations (merger driven turbulence versus 
AGN driven turbulence, and a bulk ICM in Coma versus 
condensed warm HH-gas in the NGC 1275 halo), and it is 
not fully clear in how far they are comparable. 

Shibata et al. (2001) determined the 2-point angular 
correlation function of hardness ratios as a measure of 
the temperature fluctuations detected with ASCA over an 
area of 19 square degrees in the Virgo cluster. A significant 
excess of the correlation amplitude is found at 300 kpc. 
They interpreted the random temperature fluctuations in 
Virgo-North as local heating of infalling galaxy groups. 

Future investigations should measure the pressure 
spectrum of the Coma cluster more accurately down to 
5 kpc so that the combination with the present measure- 
ments would give information about the ICM in the Coma 
cluster from 5-2800 kpc. This could give us tight con- 
straints on the type of gas turbulence, its energy content, 
the importance of magnetic fields, and on the viscosity of 
the ICM. 
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